Detecting topological order in a ground state wave function 
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A large class of topological orders can be understood and classified using the string-net condensa- 
tion picture. These topological orders can be characterized by a set of data (N,di,FliL n ,5ijk). We 
describe a way to detect this kind of topological order using only the ground state wave function. 
The method involves computing a quantity called the "topological entropy" which directly measures 
the quantum dimension D = d\ . 
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Introduction: Until recently, the only known physical 
characterizations of topological order [1] involved proper- 
ties of the Hamiltonian - e.g. quasiparticle statistics [2], 
ground state degeneracy [3, 4], and edge excitations [1]. 
In this paper, we demonstrate that topological order is 
manifest not only in these dynamical properties but also 
in the basic entanglement of the ground state wave func- 
tion. We hope that this characterization of topological 
order can be used as a theoretical tool to classify trial 
wave functions - such as resonating dimer wave functions 
[5], Gutzwiller projected states, [6-10] or quantum loop 
gas wave functions [11]. In addition, it may be useful as 
a numerical test for topological order. Finally, it demon- 
strates definitively that topological order is a property of 
a wave function, not a Hamiltonian. The classification 
of topologically ordered states is nothing but a classifi- 
cation of complex functions of thermodynamically large 
numbers of variables. 

Main Result: We focus on the (2 + 1) dimensional case 
(though the result can be generalized to any dimension) . 
Let VP be an arbitrary wave function for some two di- 
mensional lattice model. For any subset A of the lattice, 
one can compute the associated quantum entanglement 
entropy 5,4. [12] The main result of this paper is that one 
can determine the "quantum dimension" D of \t by com- 
puting the entanglement entropy Sa of particular regions 




FIG. 1: One can detect topological order in a state ^ by 
computing the von Neumann entropies Si , S2 , S3 , S4 of the 
above four regions, A\, A2, A3, A4, in the limit of R,r — * 00. 
Here the four regions are drawn in the case of the honey- 
comb lattice. The geometry ensures that the number of sites 
ni,n2,ri3,n4 along the boundaries of the 4 regions obey the 
relation m — 712 = "-3 — T14. 



A in the plane. Normal states have D = 1 while topo- 
logically ordered states have D ^ 1 (e.g. in the case of 
topological orders described by discrete gauge theories, D 
is equal to the number of elements in the gauge group) . 
Thus, D provides a way to distinguish topologically or- 
dered states from normal states. 

More specifically, consider the four regions 
Ai , A2 , ^3 , A4 drawn in Fig. 1. Let the corresponding 
entanglement entropies be Si, S2, S3, S4. Consider the 
linear combination (Si — S2) — (S3 — S4), computed in 
the limit of large, thick annuli, R, r — > 00. The main 
result of this paper is that 

{Si - S 2 ) - (S 3 - S 4 ) = - \og(D 2 ) (1) 

where D is the quantum dimension of the topological 
field theory associated with 5*. 

We call the quantity (Si - S 2 ) - (S3 - S 4 ) the "topo- 
logical entropy", — St op , since it measures the entropy 
associated with the (non-local) topological entanglement 
in "J". The above result implies that Stop is an univer- 
sal number associated with each topological phase. It is 
invariant under smooth deformations of "P. 

Physical picture: The idea behind (1) is that topo- 
logically ordered states contain nonlocal entanglement. 
Consider, for example, a spin-1/2 model with spins lo- 
cated on the links i of the honeycomb lattice and with 
a Hamiltonian realizing a Zi lattice gauge theory. [13- 
15] The ground state wave function \1/ is known exactly. 
The easiest way to describe VP is in terms of strings. One 
can think of each spin state as a string state, where a 
cr? = — 1 spin corresponds to a link occupied by a string 
and a a? = 1 spin corresponds to an empty link. In this 
language, \& is very simple: ty(X) — 1 for all string states 
X where the strings form closed loops, and ^> vanishes 
otherwise. 

All local correlations (cfcj) vanish for this state. 
However, "J contains nonlocal correlations or entangle- 
ment. To see this, imagine drawing a curve C in the plane 
(see Fig. 2). There is a nonlocal correlation between 
the spins on the links crossing this curve: (W(C)) — 
(IliGC a T) ~ 1- This correlation originates from the fact 
that the number of strings crossing the curve is always 
even. Similar correlations exist for more general states 
that contain virtual string-breaking fluctuations. In the 
general case, the nonlocal correlations can be captured 
by "fattened" string operators Wf a t (C) that act on spins 
within some distance I of C where I is the length scale 
for string breaking. 
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FIG. 2: The state \P contains nonlocal correlations originating 
from the fact that closed strings always cross a closed curve C 
an even number of times. These correlations can be measured 
by a string operator W(C) (thin blue curve). For more general 
states, a fattened string operator Wf at (C) (light blue region) 
is necessary. 




FIG. 3: A simply connected region R in the honeycomb lat- 
tice. We split the sites on the boundary links into two sites 
labeled i m and j m , where m = 1, ...,n. 



To determine whether a state is topologically ordered, 
one has to determine whether the state contains such 
nonlocal correlations or entanglement. While it is diffi- 
cult to find the explicit form of the fattened string oper- 
ators Wf a t,[16] one can establish their existence or non- 
existence using quantum information theory. The idea is 
that if the string operators exist, then the entropy of an 
annular region (such as A\ in Fig. 1) will be lower than 
one would expect based on local correlations. 

The combination [8\ — S?) — (S3 — S4) measures exactly 
this anomalous entropy. To see this, notice that (Si — 
S2) is the amount of additional entropy associated with 
closing the region A2 at the top. Similarly, (53 — S4) is 
the amount of additional entropy associated with closing 
the region A4 at the top. If 'F has only local correlations 
with correlation length £ then these two quantities are the 
same up to corrections of order 0(e~ R ^), since A2.A1 
only differ by the region at the bottom. For such states, 
limij_ >00 (S'i — S2) — (S3 — Si) = 0. Thus, a nonzero value 
for Stop signals the presence of nonlocal correlations and 
topological order. 

The universality of S top can also be understood from 
this picture. Small deformations of 'F will typically mod- 
ify the form of the string operators Wf a t and change 
their width I. However, as long as I remains finite, 
(Si — S2) — (S3 — S4) will converge to the same universal 
number when the width r of the annular region is larger 
than I. 

A simple example: Let us compute the topological en- 
tropy of the ground state wave function *F of the Z2 model 
and confirm (1) in this case. We will first compute the 
entanglement entropy Sr for an arbitrary region R. To 
make the boundary more symmetric, we split the sites on 



the boundary links into two sites (see Fig. 3). The wave 
function 'F generalizes to the new lattice in the natural 
way. The new wave function (still denoted by has the 
same entanglement entropy. 

We will decompose >F into ^ — J2i \l>} n \l/° ut where *j n 
are wave functions of spins inside R, iF° ut are wave func- 
tions of spins outside R, and I is a dummy index. A 
simple decomposition can be obtained using the string 
picture. For any qi,...,q n , with q m = 0,1, and J2 m q m 
even, we can define a wave function qn on the spins 

inside of R: vfj? q (X) = 1 if (a) the strings in X form 
closed loops and (b) X satisfies the boundary condition 
that there is a string on i m if q m = 1, and no string if 
q m = 0. Similarly, we can define a set of wave functions 
^°i*.r„ on the spins outside of R. 

If we glue * m and 1 i fOUt together - setting q m = r m for 
all m - the result is <F . Formally, this means that 

q> = V* \& in * out (2) 

9l,-*, 9l-,fa ^ > 

q 1 + ...+qi even 

It is not hard to see that the functions {'Fg" 9n : 
XL even}, and {*™. t . )I . Il : J2 m r ™ even} are or- 
thonormal up to an irrelevant normalization factor. 
Therefore, the density matrix for the region R is an 
equal weight mixture of all the 9n : Em 1 m even}. 

There are 2™ _1 such states. The entropy is therefore 
5« = (n-l)log2. [12] 

This formula applies to simply connected regions like 
the one in Fig. 3. The same argument can be applied to 
general regions R and leads to Sr = (n — j) log 2, where 
n is the number of spins along dR, and j is the number 
of disconnected boundary curves in dR, 

We are now ready to calculate the topological entropy 
associated with >F. According to (1) we need to calculate 
the entropy associated with the four regions shown in 
Fig. 1. From Sr = (n — j) log 2, we find Si = (rii — 
2) log 2, S 2 = (77 2 - 1) log 2, S 3 - (773 - 1) log 2, and S 4 = 
(774 — 2) log 2, where 771 , 772 , 773 , 774 are the number of spins 
along the boundaries of the four regions. The topological 
entropy is therefore — S% op = (ni —772 — 773 + 774 — 2) log 2. 
But the four regions are chosen such that (771 — 772) = 
(773 — 774). Thus the size dependent factor cancels out 
and we are left with — St op = —2 log 2 = — log(2 2 ). This 
is in agreement with (1) since the quantum dimension of 
Z 2 gauge theory is D = 2. 

General string-net models: To derive (1) in the gen- 
eral case, we compute the the topological entropy for the 
exactly soluble string-net models discussed in Ref. [17]. 
The ground states of these models describe a large class 
of (2 + 1) dimensional topological orders. The models 
and the associated topological orders are characterized 
by several pieces of data: (a) An integer TV - the num- 
ber of string types, (b) A completely symmetric tensor 
Sijk where i,j, k = 0, 1, N and Sijk only takes on the 
values or 1. This tensor represents the branching rules: 
three string types k are allowed to meet at a point 
if and only if S^k = 1. (c) A dual string type i* cor- 
responding to each string type i. This dual string type 
corresponds to the same string, but with the opposite 
orientation, (d) A real tensor di and a complex tensor 
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FIG. 4: A typical string-net state on the honeycomb lattice. 
The empty links correspond to spins in the i = state. The 
orientation conventions on the links are denoted by arrows. 



^kia satisfying certain algebraic relations [17]. For each 
set of Flj™ , di, Sijk satisfying these relations, there is a 
corresponding exactly soluble topologically ordered spin 
model. 

The spins in the model are located on the links k of 
the honeycomb lattice. However, the spins are not usual 
spin-1/2 spins. Each spin can be in N + 1 different states 
which we will label by i = 0, 1,...,N. Each spin state 
can be thought of as a string-net state. To do this, one 
first needs to pick an orientation for each link on the 
honeycomb lattice. When a spin is in state i, we think 
of the link as being occupied by a type-i string oriented 
in the appropriate direction. If a spin is in state i = 0, 
then we think of the link as empty. In this way spin 
states correspond to string- net states (see Fig. 4). The 
Hamiltonian of the model involves a 12 spin interaction 
[17]. The model is known to be gapped and topologically 
ordered and all the relevant quantities - ground state de- 
generacies, quasiparticle statistics, etc., can be calculated 
explicitly. 

The ground state wave function $ of the model is also 
known exactly. It is easiest to describe in terms of the 
string-net language. If a spin configuration {ik} corre- 
sponds to an invalid string-net configuration - that is, a 
string-net configurations that doesn't obey the branching 
rules defined by S^k - then &({ik}) = 0. On the other 
hand, if {ik} corresponds to a valid string-net configura- 
tion then the amplitude is in general nonzero. We would 
like to have an explicit formula for these amplitudes. Un- 
fortunately, this is not possible in general. However, we 
can write down linear relations that determine the am- 
plitudes uniquely. These relations relate the amplitudes 
of string-net configurations that only differ by small local 
transformations. The relations are given by 
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where the shaded areas represent other parts of the 
string-nets that are not changed. Also, the type-0 string 



is interpreted as the no-string (or vacuum) state. The 
first relation (3) is drawn schematically. The more precise 
statement of this rule is that any two string-net configura- 
tions on the honeycomb lattice that can be continuously 
deformed into each other have the same amplitude. In 
other words, the string-net wave function $ only depends 
on the topologies of the network of strings. 

By applying these relations multiple times, one can 
compute the amplitude for any string-net configuration 
(on the honeycomb lattice) in terms of the amplitude of 
the vacuum configuration. Thus, (3-6) completely specify 
the ground state wave function <I>. 

Let us first compute the von Neumann entropy Sr of 
the exact ground state wave function <!> for a simply con- 
nected region R (see Fig. 3). Again we split the site on 
the boundary links into two sites. We decompose $ into 
<f> = J2i $J n $° ut where $J n are wave functions of spins 
inside R, <&° ut are wave functions of spins outside R, and 
I is some dummy index. 

A wave function <E> ln on the spins inside of R can be 
defined as follows. Let {ik} be some spin configuration 
inside of R. If {ik} doesn't correspond to a valid string- 
net configuration - that is one that obeys the branching 
rules, then we define <& m ({ife}) = 0. If {ik} does corre- 
spond to a valid string-net configuration, then we define 
<!> m ({ifc}) using the same graphical rules (3-6) that we 
used for 

However, there is an additional subtlety. Recall that in 
the case of the graphical rules could be used to reduce 
any string-net configuration to the vacuum configuration. 
To fix <£, we defined $ (vacuum) = 1. 

In this case, since we are dealing with a region R with 
a boundary, string-net configurations cannot generally be 
reduced to the vacuum configuration. However, they can 
be reduced to the tree-like diagrams X^ q s y shown in Fig. 
5a. Thus, to define <& m , we need to specify the amplitude 
for all of these basic configurations. There are multiple 
ways of doing this and hence multiple possibilities for $ ln . 
Here, we will consider all the possibilities. For any label- 
ing q 1 ,...,q n ,si,...,s n - 3 of the string- net in Fig. 5(a), 
we define a wave function s j by g \{X{ q ',s'}) — 
fi{q},{q'}fi{s},{s'}- Starting from these amplitudes and us- 
ing the graphical rules (3-6) we can determine S \(X) 
for all other string-net configurations. In the same way, 
we can define wave functions ^f^ t \ on the spins outside 

of R through ®{r%( Y {r',t>}) = S {r},{r'}8{t},{t'}, where the 
Y{ r ,t} are shown in Fig. 5(b). 

Now consider the product wave functions ^V 1 ! $? ut (1 . 
These are wave functions on the all the spins in the sys- 
tem - both inside and outside R. They can be gener- 
ated from the amplitudes for the string-net configurations 
z { q ,s,r,t} in Fig. 6: 

®Tq,s}®{rl}( Z W,s>,r',t>}) = ${l};W} S {s},{s'} S {r},{r'} S {t}.,{t'} 

On the other hand, it is not hard to show that for the 
ground state wave function $, the amplitude for Z^ q s r ^ 
is 

®{Z{ q ,s,r,t}) = S {Q},{r} S {s},{t} HCV^sT) 
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FIG. 5: The basic string-net configurations (a) Xi qiS \ for in- 
side R and (b) Ys r j\ for outside R. 




FIG. 6: The string-net configuration ^{ q , s , r ,t} obtained by 
"gluing" the configuration X{ q ^ s y to the configuration Kr ri t> 
in Fig. 5. 

Comparing the two, we see that 

{q,s,r,t}, m 

It turns out that the wave functions {^/^ a i} are or ~ 
thonormal, as are the {${"*■}.} (up to an irrelevant nor- 
malization constant). This means that we can use them 
as a basis. If we denote $™ q s } by \{q, s, r, t}), then 
in this basis, the wave function <E> is 

({q,s,r,t}\$) = ${ q },{r}S{ s },{t}Y[(VdZ) ( 8 ) 

m 

The density matrix for the region R can now be obtained 
by tracing out the spins outside of R, or equivalently, 



tracing out the spin states \{r, t}): 

(W, s'}\pR\{q, s}) = <*{„},{«'}<*{«},{*'} JJ d qm (9) 

m 

Since the density matrix is diagonal, we can easily obtain 
the entanglement entropy for Sr. Normalizing pn so that 
Tr(pu) — 1, and taking — Trpn log pr, we find 

where D = d\. The sum can be evaluated explicitly 
(with the help of the relations in [17]). The result is 

5 fl =-log( J D)-n£§log(|) (11) 

This result applies to simply connected regions like the 
one shown in Fig. 1. The same argument can be applied 
to general regions R. In the general case, we find 

AT fA\ 

5* = -,-tog(D)-n£%l°g(%) (12) 

where n is the number of spins along dR, and j is the 
number of disconnected boundary curves in dR. 

We can now calculate the topological entropy associ- 
ated with Applying (12), we find Si = -2\ogD - 
TH so, S2 = -log-D - n 2 so, S3 = -logD — n 3 so, and 
54 = —2 log D — 7J4S0 where ni, 712,^3,714 are the num- 
bers of spins along the boundaries of the four regions, and 

so = J2k=o ~d {%)• ^ n e topological entropy is there- 
fore —Stop — -2 log£) + (ni — n 2 — n 3 + n 4 )s = — 21ogD 
in agreement with (1). 

Near the completion of this paper, we become aware of 
a similar result, obtained independently in the recent pa- 
per, Ref. [18]. This research is supported by NSF Grant 
No. DMR-04-33632 and by ARO Grant No. W911NF- 
05-1-0474. 
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